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In Ref. 1 we have mistakenly claimed that the applicability of the Hellmann-Feynman theorem in fixed-node quantum 
Monte Carlo calculations is not subject to-.the manner how the nodal boundary depends on an external parameter A. 
As it has been pointed out by Huang et al.0 in their comment on Ref. 1, this statement is not correct in general, except 
1 1 , where the Hellmann-Feynman force is calculated for a nodal boundary which coincides with that of the unconstrained 
exact eigenfunction. We want to point out the error in our arguments and present an explicit expression for the 
correction term which supplements the Hellmann-Feynman force. _ 

In our approach the fixed-node approximation is treated as a Dirichlct type of boundary value problemtl on a nodal 
region Q. Properly stated O has to be taken as an open subset of R 3N in which the fixed-node wavefunction ^(A, fi) 
g . satisfies the Schrodinger equation 

#(A)*(A,n) = £?[A,n]*(A,n). (i) 

rH , 

The fixed-node wavefunction 'J (A, f2) has to vanish on the boundary <9f2 of f2 and must be continous on the entire 
q ' 57 U dCl. In order to define the derivatives of the wavefunction on the boundary dil one has to take the limit of 
derivatives of interior points when approaching the boundaryu. Within such a framework no 5 function term appears 
for the second derivatives on dfl since we do not extend our wavefunction beyond the boundary. Neclecting spurious 
\ singularities in the potential V(A) on dil we can actually conclude from 

' \ 

" A* A,Q =2 F A -BA,Q )*(A,fi) (2) 

(N : v ' 

that the limits of the second derivatives vanish almost everywhere on the boundary. We want to stress however that 
this is not in contradiction to the arguments given in Ref. 2, where the wavefunction has been extended over the 
whole space. In this case one actually encounters discontinous first derivatives when crossing the nodes. 

Revisiting Eq. (8) of Ref. 1, we can identify the missing boundary term mentioned in Ref. 2 by inspection of the 
second line. After differentiation with respect to A we obtain 



/ [0Atf(O)flo*(O)+*(O)0 A flA *(0) + *(0) H 8x9(0)] cM (3) 
Jn(o) 1 A=0 J 

... n 

The first term in the integral vanishes due to Eq. ( |1|) and the normalization constraint Eq. (6) in Ref. 1 from which 

O ■ follows E 



, *(Q) d x 9(0) dn = 0, (4) 
U ' 

whereas the second term yields the standard Hellmann-Feynman force, however the third term does not vanish in 
general as we have erroneously claimed in Ref. 1. Applying Green's second formulaa we can rewrite the third term 
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tf (0) H 3 A *(0) dQ, = -- / I V*(0) I 9 A *(0) ddil (5) 
n(o) 2 J9Q( ) 



x To see that this is valid for a parameter dependent boundary we refer to Ref. 1 for a discussion of the boundary terms. 
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where the sign of the boundary term corresponds to \& > inside the nodal region. For this step we presume that 
the first and second derivatives of ^(0) and d\ ^(0) can be continuously extended to the boundary in the sense 
discussed above. In order to get a better understanding of the physical character of this term we have to generalize 
our considerations by allowing the external parameter A and the nodal domain to vary independently. Doing so we 
can rewrite the derivative 8x^(0) like a total differential 

9a*(0) = [ d A *(A, O(0)) + 9 A *(0, 0(A)) ]| A=0 (6) 

The first term ^^(A, O(0)) corresponds to the change of the wave function with respect to the external parameter A 
under the constraint that the nodes are kept fixed. It vanishes on the boundary and therefore does not contribute to 
Eq. ( ||). The second term represents the change of the wave function ^^(0, 0(A)) under a variation of the nodes 
only. Since ^(0, 0(A)) is an eigenfunction of Ho on the nodal domain O(A) we obtain 

Ho d A *(0, O(A))| a=0 = d x E[0, 0(A)]| A=0 *(0, 0(0)) + E[Q, 0(0)] A *(O, O(A))| A=0 (7) 

where the contribution of the second term to the left side of Eq. ( |^) vanishes due to the normalization contraint 
( ^). The modified Hellmann-Feynman theorem for fixed-node quantum Monte Carlo calculations with parameter 
dependent nodal boundary is therefore given by 

9a£[A,O(A)]| a=0 = / *(0)dxH x *(0)rfO (8) 
Jnio) x=0 



0(0) 

|w(o)| a A *(o,o(A))| A=0 ddn 

3O(0) 

*(0) d\H\ (0) dn + d x E[0, 0(A)] | A=0 
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with an additional term which can be interpreted as the linear response of the energy with respect to the variations 
of the the nodal region. This term vanishestJ provided that the nodal region 0(0) coincides with a nodal region of the 
unconstrained solution of the Schrodinger equation. 

Finally we want to mention that the main part of our paper concerning the dipole moment of CO is not affected 
by this correction. In our actual calculations we have only used parameter independent nodal boundaries for which 
the unmodified Hellmann-Feynman theorem remains applicable. However, further studies of the implications of the 
additional term seem to be necessary. 
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